Universality in escape from a modulated potential well 
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We show that the rate of activated escape W from a periodically modulated potential displays 
scaling behavior versus modulation amplitude A. For adiabatic modulation of an optically trapped 
Brownian particle, measurements yield In oc (Ac — A)'' with fi — 1.5. The theory gives fi = 3/2 in 
the adiabatic limit and predicts a crossover to /i = 2 scaling as A approaches the bifurcation point 
where the metastable state disappears. 
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Activated processes, such as escape from a metastable 
state, are exponentially sensitive to time-dependent ex- 
ternal perturbations. Escape in driven systems has been 
investigated in studies of Josephson junctions ||l|, |j and 
infrared photoemission |3| . It has been also explored re- 
cently in the context of thermally activated magnetiza- 
tion reversal in nanomagnets driven by a time-dependent 
magnetic field or a spin-polarized current The 
strong dependence of the rates of activated processes on 
the wave form of the driving field enables selective control 
of the rates, which is important for many applications, 
for example control of rate and direction of diffusion in 
spatially periodic systems 

In this paper we investigate the dependence of the rate 
of activated escape W on the amplitude A of a periodic 
driving field. We show that this dependence displays 
universal behavior for large A and investigate the scal- 
ing of In W with A. Experimental results are obtained 
for a colloidal particle trapped in a modulated optically- 
generated potential. This system was used earlier Q for 
a quantitative test of the Kramers theory Q of activated 
escape over a stationary potential barrier. We observe 
a power-law behavior of In 14^ as A approaches a critical 
value. The critical exponent for adiabatic modulation is 
fi = 1.5, in agreement with theory. Adiabaticity is always 
violated close to the bifurcation point, and we predict the 
emergence of a different scaling, with fi — 2. 

Field dependence of the escape rate is well understood 
for adiabatically slow modulation, where the field fre- 
quency ujp is small compared to the reciprocal relax- 
ation time . The driven system remains instanta- 
neously in thermal equilibrium. The adiabatic escape 
rate Wad oc exp[— i?ad/fcsr] depends on the field primar- 
ily through the instantaneous activation energy i?ad(^)- 
Even where the modulation of i?ad is small compared to 
the zero-field i?ad, it may still substantially exceed fcsT, 
leading to strong modulation of VFad- 

Finding the escape rate becomes more complicated for 
nonadiabatic modulation, where u/p ^ t^^^ because a 
driven system is far from thermal equilibrium. The prob- 
lem has been addressed theoretically for different models 
of fluctuating systems §, |, |l^, |l|, |l| . The underly- 
ing idea is that activated escape results from an optimal 
fluctuation, which brings the system from a metastable 
to an appropriate saddle-type periodic state. We use it 
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FIG. 1: Amplitude scaling of the period-averaged rate of es- 
cape W from the modulated potential well shown in the in- 
set {Wo is the escape rate in the absence of modulation). 
The metastable state of forced vibrations disappears once the 
modulation amplitude reaches the critical value Ac- 



here to draw general conclusions about the scaling of the 
period-averaged escape rate W with the field amplitude 
A. 

The field dependence of \nW for W ^ t^^,ujp is 
sketched in Fig. 0. For small A, the field is a pertur- 
bation. Its major effect is to heat the system, with 
the change of the effective energy-dependent tempera- 
ture quadratic in A (the lowest order term after period 
averaging). As a result, In [W/IVo] oc A^, where Wq is 
the escape rate in the absence of modulation |0, ||. The 
range of A is limited to In [W/Wq] ^ 1. 

For stronger modulating fields, the activation energy 
of escape i? oc In W becomes linear in the amplitude 
A jl^. This behavior occurs because the optimal fluc- 
tuation leading to escape is a real-time instanton [p^ , 
with duration ^ t^. In the absence of modulation, es- 
cape has the same probability to happen at any time. 
Modulation lifts the time degeneracy and synchronizes 
escape events. As a result, even for zero-mean modula- 
tion In W/H/q oc a. For an overdamped Brownian parti- 
cle, an explicit solution was obtained ||ll| throughout the 
heating and log-linear regions in Fig.^ and the results 
have been confirmed by extensive simulations [T^ . 

As we show in this paper, there is another region where 
the dependence of W on A is universal and unexpected: 
the vicinity of the saddle-node bifurcation point Ac where 
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the metastable state of forced vibrations merges with a 
saddle-type periodic state [|6|. In this range one of the 
motions of the system becomes slow. The universality is 
related to the corresponding critical slowing down. How- 
ever, it turns out that the system can display different 
types of critical behavior depending on the relationship 
between ujp and the relaxation time in the absence of 
driving t^o- 

Experiments on amplitude-dependent escape from a si- 
nusoidally modulated potential were performed using an 
optically-trapped Brownian particle in water. A config- 
urable optical potential was constructed by focusing two 
parallel laser beams through a single microscope objec- 
tive lens. Each beam creates a stable three-dimensional 
trap as a result of electric field gradient forces exerted on 
a transparent dielectric spherical silica particle of diam- 
eter 0.6 /zm. At the focal plane, the beams are typically 
displaced by 0.25 to 0.45 //m, creating a double-well po- 
tential. The depths of the potential wells and the height 
of the intervening barrier, which is typically 1 to 10 fc^T, 
are readily controlled by adjusting beam intensities and 
separations. 

The two HeNe lasers (17 mW, 633 nm) that create the 
traps are stabilized by electro-optic modulators and im- 
aged into a sample cell by a lOOx objective. The beams 
are mutually incoherent and circularly polarized when 
they enter the microscope. A single trapped colloidal 
sphere is imaged onto a digital camera operating at 200 
frames/s. The x and y coordinates of the sphere's cen- 
ter are computed using a pattern matching routine that 
yields spatial resolution better than 10 nm, whereas the 
z coordinate parallel to the light wave vector is extracted 
by analysis of the image as the particle position fluctuates 
about the focal plane. The overall analysis and storing 
of the sphere's (x, y, z) coordinates is accomplished in 
less than a frame duration so that no images need to be 
recorded. 

The stability of the system is sufficient to compute the 
full 3-dimensional optical potential from long-time mea- 
surements of the particle probability density 0|. The 
absolute transition rates calculated from the Kramers 
expression using measured curvatures at the stable and 
unstable points of the potential are in quantitative agree- 
ment with the experimental data over 4 orders of magni- 
tude. 

The electro-optic modulators that stabilize the optical 
traps also enable ac modulation with arbitrary waveforms 
derived from an electronic signal generator. When the 
two beams overlap at the separations required to form 
a potential barrier in the range 5 to 8 fcsT, each con- 
tributes to the position and shape of both potential min- 
ima and the potential in the intermediate range. As a 
consequence of this nonlocality, a small modulation of 
one beam causes a nearly antisymmetric change of both 
barrier heights. 

The results on over-barrier transitions from one of the 
wells for the modulating signal F{t) = Aco&ujpt are plot- 
ted in Fig. ^. The amplitude A in Fig. ^ is normalized 
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FIG. 2: Modulation amplitude dependence of the maximal 
transition rate Wm normalized to Wo, the unmodulated tran- 
sition rate. The data points represent the transition rates for 
an unmodulated barrier AUo/ksT = 5.8 {Wo = 0.069 s"^). 
The solid line is a least-squares fit to the data with an asymp- 
totic form described by the theory (see text). The modulation 
frequency ujf/2tt = 0.4 Hz. Inset: the number of transitions 
for A = 2.75fc_Br versus time. 



to fcsT, where T = 298 K. The intrawell relaxation rate 
in the absence of modulation t~Q « 10'^ s~^ is much 
larger than uip, so the data fall into the adiabatic regime. 
In this case escape is most likely to occur when the in- 
stantaneous potential barrier AU{t) is minimal. For large 
modulation amplitude and for temperatures used here, 
the maximal escape rate Wm becomes comparable to the 
modulation frequency, and therefore it is Wm rather than 
the period-averaged rate W, that is the more relevant 
measure. In obtaining Wm we used the distribution of 
the dwell time illustrated in the inset to Fig. |2[ The 
dwell time was computed for the modulated double-well 
system as the time interval between the appearance in, 
and escape from, a well. 

The solid line in Fig. ^ represents a least-squares fit to 
the data for large A using the function In Wm/Wo — Ci — 
C2{Ac — A)^^, where /i = 1.5, and the critical modulation 
amplitude Ac and Ci_2 are fitting parameters fl^ ]. 

To interpret these results we will consider activated 
escape in the system for A close to the critical amplitude 
Ac where the metastable and unstable periodic states 
merge. In the adiabatic approximation, the stable and 
unstable states correspond to the potential well and the 
barrier top. Once per period, when the driving |-F(t)| is 
a maximum (for example when t = nrp, n — 0, ±1, . . .), 
the stable and unstable states are closest to each other 
and the barrier height AU(t) is a minimum. We take 
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F(0) positive, F{Q) = A. The well and barrier top merge 
for A — Ac, which is the adiabatic bifurcation point. 

For small 6 A = A — Ac and \t — nTp\, the dynamics of 
the system near the metastable state is slow [|l6[ . It is de- 
scribed by the "soft mode" with coordinate q{t). For neg- 
ative SA the system has a stable and an unstable state, qa 
and Qb- We set g = at the saddle-node state where they 
merge in the limit of small u p for A = Ac,t = nrp . Close 
to Ac and nTp the equation of motion for g, neglecting 
fluctuations, takes the form 



dq 



q^ + Fit) - Ac 



(1) 



(tro is the relaxation time for F = 0, ujptj-Q ^ 1). In the 
adiabatic approximation, F(t) depends parametrically on 
t. Activated transitions occur over the instantaneous bar- 
rier AU{t) = (4/3) [Ae - F(t)]3/2. 

The relaxation time of the system — iro[^c ~ 
F(t)]~^/2 becomes large close to the bifurcation point. 
Therefore even a slowly varying field will ultimately be- 
come too fast for the system to follow adiabatically. To 
allow for nonadiabaticity we expand F(t) in the region of 
small u!p\t — nTp\. The equation of motion (|l|) in scaled 
variables Q ^ q/il, t ~ fl{t — nTp)/tro is then 



Q^K{Q,t), K = Q^ + C-t^ 



(2) 



with C n-^iA - Ac). In Eq. @ Q = dQ/dr, and 
n = (— it^pd^F/di^)^/^ (the derivative is evaluated at 
t = and A = Ac). For sinusoidal modulation — 

The dynamics is determined by one dimensionless 
parameter C. For — C ^ 1, the adiabatic approximation 
applies, and the stable and unstable states are given by 
g„ b(r) = =f(-C + T-Y^^- On the other hand, for |C| < 
1 the time dependence of the states becomes distorted 
and asymmetric, cf. inset in Fig. ^. For C = 1, the 
states given by Eq. ^ merge and Qa{T) — Qb{'r) — 
T. The corresponding nonadiabatic bifurcational value 
of the amplitude is Ah ^ Ac + ^^ . 

Escape from the attractor can be quite generally de- 
scribed by adding a random Gaussian force to Eq. (|^). 
Because the system dynamics is slow, this force is effec- 
tively (5-correlated (white noise), with intensity D <x ksT 
for thermal noise. The corresponding force ^(r) in Eq. (||) 
for Q has a correlator (C(t)C(t')) = 2DS{t ~ r'), D = 

D/n\o.^ 

The noise-driven system is most likely to escape during 
the time when F{t) is close to its maximum. The total 
probability to escape over this time is Wt = Wtp, where 
W is the period-averaged escape rate. We will calculate 
Wt assuming that the noise intensity D is the smallest 
parameter of the theory, in which case Wt oc exp{—R/D). 
The calculation can be done by generalizing the instanton 
technique ||l5| to the Langevin equation Q = K{Q, t) + 
^(r). The activation energy of escape R = {D/D)S is 
given by the solution of the variational problem 
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FIG. 3: The scaled activation energy of escape 5 « —D In Wt 
(dots), obtained numerically from the variational problem 
(^), vs. — Q~^{A — Ac). In the adiabatic approximation 
S = |(— C)'^''^ (the dashed line), which explains the results 
in Fig. H. The nonadiabatic bifurcation point is = 1. The 
solid lines show the corrected adiabatic expression for large- 
ICIi Eq- (§)i fl'iid the nonadiabatic behavior S oc (1 — i^)^, 
Eq. (H). Inset: the attracting and repelling states Qa and 
Qb (thin lines) for weakly nonadiabatic driving, Eq. (□) with 
— 0.95, and the most probable escape path (thick line). 



The boundary conditions are Q{t) Qa{T) for t ~oo 
and Q{t) Qbi^) for r — > oo. They correspond to the 
picture in which, prior to escape, the system is perform- 
ing small fluctuations about the metastable state Qa- Es- 
cape occurs as a result of a large fluctuation that brings 
the system to the saddle state Qb- The most probable 
among such fluctuations is the one in which the sys- 
tem moves along a trajectory [the most probable escape 
path (MPEP)] that minimizes the functional (^. From 
Eq. (0), the activation energy i? is a function of one pa- 
rameter C, shown in Fig. ^ 

An explicit expression for R can be obtained in the 
adiabatic limit — C 3> 1. Here, one can ignore the term 

in K (H) for typical |t 
is an instanton Qopt{T,TQ) 
centered at an arbitrary tq. The adiabatic activation 
energy is i? oc 5 = (4/3) (-C)^/^ oc [Ac-A)^/^ ^ This 
agrees with the experimental data shown in FigT^. 

The time-dependent term in K lifts translational in- 
variance of the instanton Qopt(''", tq) and synchronizes 
escape events. The action S is minimal for tq — 0. To 
lowest order in 1/C 



< C~^/^. Then the MPEP 
|Cr/2tanh[|C|i/2(r-ro)] 



5«(4/3)|Cp/2 + (^V6)ICr'/'. 



(4) 



dT{Q - Kf 



(3) 



From (^, the nonadiabatic correction to the activation 
energy diverges as |(5A|~^/2 as A approaches Ac. 

The activation energy may also be found explicitly 
close to the bifurcation point C, = 1. In this case escape 
is most likely to occur while the coexisting attracting 
and repelhng trajectories Qa{T) and Qb{T) stay close to 
each other. Because of the special structure of the states 
Qa,b{T), the variational problem (0) for the MPEP can 
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be linearized and solved, giving 

5=(V8)i/2(l-Cf, 1-C«1- (5) 

Eq. (j^) shows that, near the nonadiabatic bifurcation 
point Ay,, there emerges another scaling region, where 
the activation energy R oc (Ab — A)"^. The nonadiabatic 
exponent is thus n = 2, in contrast to = 3/2 in the 
adiabatic limit. It describes R{A) for the modulation 
amplitude lying between Ac and ^b, as seen from Fig. |^. 
The width of this interval is cx o;;^ provided LUpty-Q <C 1. 

The above theory does not apply in the exponentially 
narrow range 1 — ^ ^ exp[—C /ajptj-o] (C ~ 1). This is 
because, as A approaches A^,, the relaxation time of the 
system becomes logarithmically long. Then the expan- 
sion of the driving force int/Tp that leads to the equation 
of motion (||) becomes inapplicable. For uiptr ^ 1 the 
duration of the MPEP greatly exceeds the period Tp . In 
this case the activation energy again displays the fi = 3/2 
power-law dependence on the distance to the true bifur- 
cation value of the amplitude The — A)'^/^-law 
also applies asymptotically if the driving is nonadiabatic 
for small A. 

The experimental data discussed above refer to very 
small ujpt^o < 10^^. Therefore the nonadiabatic region 
of modulating amplitudes is narrow and the activation 
energy there is less than ksT. Higher modulation fre- 



quencies will be required to detect the crossover from 
the observed fi — i/2 scaling to the nonadiabatic fi = 2 
scaling. It appears that the ii — 3/2 scaling function de- 
scribes the data in a broad range of A. It thus provides 
a good interpolation of the activation energy from the 
linear in A to the critical region. 

In conclusion, we have identified several regions where 
the activation energy R of escape from a metastable 
potential displays scaling behavior as a function of the 
amplitude A of an externally applied periodic modula- 
tion. We have demonstrated with a system whose poten- 
tial was directly measured that, near the amplitude Ac 
where the local minimum and maximum of the potential 
contact one another, R cx {Ac — A)'^/^. We have also 
shown that, because of nonadiabaticity near the bifurca- 
tion point, there necessarily emerges a region where the 
scaling exponent changes from n = 3/2 to fi = 2. We ex- 
pect that these scalings can be observed in other systems, 
such as modulated Josephson junctions and nanomag- 
nets 1^, ||, |l9|. The ideas discussed here can be used to 
gain additional insight into the physics of these systems 
and more generally, activated effects in nonequilibrium 
systems. 
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